What happens to geometric phase when spin-orbit interactions lift band degeneracy? 
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Spin-orbit interaction lifts accidental band degeneracy. The geometric phase 7(C) — ±7r for 
circuits surrounding a line of such degeneracy cannot survive completely unchanged. The change 
depends on how the spin is fixed during adiabatic evolution. For spin fixed along the internal spin- 
orbit field, 7(C) decreases to zero as the circuit collapses around the line of lifted degeneracy. For 
spin fixed along a perpendicular axis, the conical intersection persists and 7(C) = ±7r is unchanged. 

PACS numbers: 71.70.Ej, 71.18.+y 



Geometric (or Berry) phased has become a power- 
ful tool for analysis of waves in periodic systems, espe- 
cially electrons in crystals^^^. Wavevector k pro- 
vides a space in which adiabatic evolution of wavefunc- 
tions ipn(k, r) can be studied. Singular behavior occurs at 
band degeneracies where energies ei(fe) = 62(h) are equal. 
In crystals with inversion symmetry, ignoring spin-orbit 
interactions, degeneracies occur along closed lines in k- 
space^. The periodic part u n (k, r) — exp(— ik ■ r)ijj n (k, r) 
of tp is an eigenstatc of TL(k) = (p + hk) 2 /2m + V(r). 
Let the wavevector k(t) be given a time evolution which 
takes it on the circuit C, with k(T) = k(Q). Now 
suppose that wavefunction evolution is determined by 
the time-dependent Schrodinger equation with the time- 
dependent Hamiltonian H(k(t)). The time-evolution 
is assumed adiabatic, namely u n (k,r,t) cx u n (k(t),f). 
Berry's argument shows that u n (k,f,T) differs from 
u n (k, f, 0) by the factor exp[«7(C, T)], where the phase 
7(C, T) has two parts, 7(C) + j(T). The dynamical part 

j(T) = — L dte n (k(t))/Ti depends on the time elapsed, 
and the geometric part 
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is invariant and intrinsic to the circuit and the band 
properties. In particular, 7(C) = ±tt if C encloses one 
(or an odd number) of degeneracy lines. This change of 
wavefunction sign is familiar from other problems where 
a circuit of adiabatic evolution surrounds a conical in- 
tersection. Direct evaluation of Eq.([T]) is problematic. 
Wavefunctions must be continuous and single-valued. 

Although gauge invariance is not evident in Eq.([TJ), 
Berry gave also an alternate form, for a 3-dimcnsional pa- 
rameter space fc, as the flux through a surface S (bounded 
by C) of a vector V n . 



7 (C) = - / dS s • V n 
Js 



V n = ImJ2 



(n\V R H\m) x (m\V s H\n) 
[e m (k) - e n (k)} 2 
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The gauge invariance of this vector is evident. Conditions 
of continuity and single- valuedness of wavefunctions are 
no longer required. If the circuit surrounds a singular- 
ity described by a 2 x 2 effective Hamiltonian, then the 
flux equals half the solid angle f2(C) subtended in an 
appropriate scaled space by the circuit as seen from the 
point of singularity. The appropriate scaled space is the 
one in which the 2x2 Hamiltonian for states near the 
conical intersection has the form Ti c ff = R ■ a in terms 
of scaled coordinates R — (X, Y, Z) and Pauli matrices 
<? = (&x, VyiVz)- This method will be used twice in this 
paper. The eigenvalues are ±i? = p z R, where the quan- 
tum number p z — ±1 is introduced as a branch index. 
The geometric phase is then 7(C) = —p z il(C)/2. 

Mikitik and Sharlai 7 - provide convincing evidence that 
the geometric phase ±7r is seen experimentally as a shift 
in the semiclassical quantization condition^ determin- 
ing the de Haas- van Alphen oscillations. An extreme 
experimental case is the shifted quantum Hall oscilla- 
tions originating from orbits near the "Dirac points" in 
graphen o 11 ' 12 . The shifts of quantization condition occur 
for electron orbits (in a i?-field) which surround a degen- 
eracy line (or point, for graphene.) They also argue 10 
that spin-orbit effects can mostly be ignored. This 
is correct for lighter elements with spin-orbit strength 
^/A<1, A being any other relevant electron scale such 
as a band gap. However, the mathematics and the cor- 
rections need elucidation. Spin-orbit coupling destroys 
band degeneracy lines. It is not evident what happens to 
the geometric phase of ±7r. 

To see the effect of spin-orbit interactions, add to TC(k) 
the piece Hso = (<?/4m 2 c 2 ) • VV^ x (p + hk). Choose 
some point k* of accidental degeneracy, and find ener- 
gies and eigenstates at nearby fc-points using degener- 
ate k ■ p perturbation theory. For notational simplicity, 
k* is the temporary origin of k. The degenerate basis 
functions |1) and |2) are the periodic parts u\ and ui at 
k = k*. A phase convention is needed; the coefficients 
Cq of the expansion u(r) = Cg exp(iG • f) are chosen 
real. This requires inversion symmetry, which is here- 
after assumed. Each state has two spin orientations, so 



the effective Hamiltonian matrix is 4 x 4, with the form 

^ ^ _ ( hk-v a l hk-Vbl-i^-a\ 
\ hk ■ Wftl + i£ ■ a —hk ■ v a l J 

where 1 and t? are 2x2 matrices in spin space. Terms 
proportional to the 4x4 unit matrix do not mix or split 
the states and are omitted. The vector v a is (vi — W2)/2, 
where v n is the band velocity Vfce n /?i at the degener- 
acy k*. The vector Vb is the off-diagonal term (2|p/m|l), 
which is pure real since Cq is real. The vector it; is 
(2|Vy x p|l)/4m 2 c 2 . This is pure imaginary since there 
is also time-reversal symmetry, under an assumption of 
no magnetic order or external B-field. Thus three real 
vectors, v a , Vb, and £, determine the bands near k* . The 
vector £ is a close analog to angular momentum. Con- 
sider a system with two degenerate p-states \x) and \y). 
The angular momentum operator L has an imaginary off- 
diagonal element. The mixed states |a;) ± i\y) are eigen- 
states of L with (L) = ±mhz. The magnitude m deviates 
from 1 if the point symmetry is less than spherical. The 
vector £ will be called the "orbit moment." 

First suppose that £ = 0. Since v a and Vf, are not 
generally co-linear, they define a direction of k, namely 
v a X Vb, along which Hcs = 0. This is the direction of 
the line of degeneracy. After allowing £ ^ 0, eigenvalues 
of Eq.Q are ±A where 

a = ^i + ^i+e (5) 

with K a = Hk ■ v a , Kb = hk ■ Vb, and £ = |£|. Each eigen- 
value belongs to a Kramers doublet of two opposite spin 
states. The original degeneracy (without spin-orbit inter- 
action) of 2 (neglecting spin) or 4 (including spin) is lifted 
everywhere unless £ = 0. This should happen only at 
isolated points in the Brillouin zone, not coinciding with 
degeneracy lines k*. No accidental degeneracies remain, 
but Kramers degeneracy occurs everywhere. Bands near 
k* are shown in Fig[TJ 

The geometric phase under consideration involves 
a circuit C(k) surrounding the k* line. A circular 
path in two-dimensional (re a , Kf,)-space, namely C — 
(ncoscj), Ksin0), — > <f> — > 2n is the simplest realiza- 
tion. To calculate 7(C), separate Eq.(j4]) into two similar 
2x2 submatrices by choosing basis states with spins po- 
larized along £, which will be used as the z-axis of spin 
space. The submatrices are 

H «-\K h ±it -K a J' (6) 

where the upper sign goes with spin up, a z = 1. 

The circuit can now be considered as a path C(A) in 
a 3-d A-space, where (X x , X y , \ z ) = (Kb, <7j£, K a ). On this 
circuit, A, k, and £ are all constant. The effective Hamil- 
tonian has the desired scaled form. The solid angle is 




FIG. 1: Energy versus \k\ near the degeneracy point, for (a) 
no magnetic field, (b) field parallel to £, and (c) field perpen- 
dicular to f. The dashed lines are for £ = and b — 0; solid 
lines in panel (a) are ±A, which becomes ±£ at the degeneracy 
point k — 0. 

<7 z 27r(l — £/A), so the geometric phase is 

7 (C) = -(A z ^)7r(l-£/A). (7) 

where A z = ±1 is the branch index. This is one of the 
two main results of this paper. It shows how spin-orbit 
splitting destroys the simple phase of ±7r when the cir- 
cuit has such a small radius that £ ~ A. If spin-orbit 
interaction is weak, it does not need a large orbit to have 
£/A <C 1 and approach the full simple phase of ±tt. 

This is not the full story. The choice to evolve at fixed 
a z was arbitrary. The states of Kramers doublets can 
be mixed by arbitrary unitary transformations. Evolu- 
tion of a doublet around a circuit introduces not a sim- 
ple geometric phase, but a unitary matrix. The 7(C) 
phases just computed are actually the diagonal elements 
exp(±i7(C)) of a 2 x 2 unitary matrix in the represen- 
tation with spin quantized along £. It will emerge below 
that this is indeed the correct adiabatic evolution of the 
Kramers doublet when an small magnetic field is imposed 
along the £ direction. 

Berry's original argument assumed that TL had a dis- 
crete spectrum along C. There is a physically natural 
way to retain this. Magnetic fields present in experiment 
since they are used to cause cyclic evolution in fc-space. 
Magnetic fields also lift Kramers degeneracy. The sim- 
plest theoretical device is to add to TL c s a Zeeman term 
7iz = — b ■ a coupling only to spin. 

To proceed further, an explicit representation of eigen- 
states is needed. Eigenstates of the effective Hamiltonian 
©, labeled by energy ±A and a z =1,1 are chosen as 

l-> = l-M> = i(-^+ A <e )®IT>, (8) 
l«> = l + M> = i(£+£)®lt>, do) 
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i«>=i+m>=£(£:^)®u>. (id 

These are written as direct product of spatial times spin 
two-vectors. The normalization is n = y2A(A + n a ). As 
long as £ is non-zero, 1/n is non-singular and these are 
smooth, single- valued functions of (re a , Kb), unique except 
for an arbitrary overall phase, which cannot alter 7(C). 
The lower Kramers doublet |s), \t) has "orbit moments" 
(i|VV x p/4m 2 c 2 |i) = T(£/A)£ oriented antiparallel to 
spin, while the upper Kramers doublet \u), \ v) has iden- 
tical orbit moments except oriented parallel to spin. 

Now the Zeeman term is added. Diamagnetic coupling 
is neglected. Without loss of generality, the part of the 
field B — b/fiB perpendicular to £ can be used to define 
the x direction of spin. The total Hamiltonian in the 
basis \s), \t), \u), \v) is 

/ X + b z fb x e iu i^b x \ 
v = _ fb x e-^ X-b z -ifa 
tot ifa -X + b z fb x e-™ 

V -i$b m fb x e^ -X-b z ) 

(12) 

The factor (re/A) exp(za>) = (s|<r+|t) introduces the new 
angle u> 

#> = ~ - ^P\- (13) 

K K(A + K a ) 

As the circuit C is followed (</> going from to 2ir, with 
£, k, A constant), u) also evolves from to 2ir. 

If the field b is along z, the upper and lower Kramers 
doublets are not coupled. The degeneracy is lifted ev- 
erywhere, and adiabatic evolution proceeds smoothly on 
the resulting non-degenerate states, yielding the phases 
7(C) ofEq.([7]). The previous discussion was correct. The 
result [7] can also be obtained directly from Eq.([T]) using 
Eqs. (|8l9i UOIlip . For fields perpendicular to z, there is 
both intra- and inter-doublet spin mixing, according to 
Eq. (fT2")l . To first order, since b <C A, inter-doublet mix- 
ing terms ±££6 X /A can be neglected, giving 2x2 effective 
Hamiltonian matrices, of the form 

n^)=x z xi-[ fb J^ ^_ e bz ) (i4) 



The eigenvalues are 

K 2 

±\±fj, where ^ 2 = b 2 z + -^b 2 x (15) 

A 

These eigenvalues have an interesting feature: at the de- 
generacy point k = 0, in the center of circuit C, \x = and 
Kramers degeneracy is not lifted, provided b is perpen- 
dicular to £. The states at k* have anisotropic g factors 
which vanish in two directions. The vanishing Zeeman 
splitting means that a conical intersection, hidden unless 
6/0, exists exactly where the original band intersec- 
tion (for £ = 0) was located. This also yields a simple 
geometrical phase of ±7T. Bands for b \\ £ and b _L £ are 
shown in FigQ] panels (b) and (c) . 

A full calculation of 7(C) for the 4 new eigenstates 
of Eq. (fl"2"|) is difficult. The Berry method of solid an- 
gle works when the basis functions |1), \2) of the 2x2 
effective Hamiltonian are fixed at k*, whereas the basis 
functions \s), \t) or \u), \v) used in Eq. (fT4"]) depend on k. 
However, the most important limit remaining to be re- 
solved is when the circuit radius k is small relative to 
spin-orbit splitting £. In this limit, the basis functions 
loose their fc-dependence. The relevant scaled parame- 
ters are jl = ((n/X)b x cosw, — A 2 (k/A)& x sinw, b z ). The 
circuit parameterized by <j) is equally well parameterized 
by uj which evolves from to 2-7T. The solid angle in 
/2-space is 27rA 2 (l — b z /fJ>), so the geometric phase is 

7 (C) = -7rf3 z A z (l - -^=J^=] , (16) 

where (A z ,j3 z ) are the two branch indices in the eigen- 
value ±A±/i = A z A+/3j/i. This is the other main result of 
this paper. If b z = 0, the full geometric phase 7(C) = ±ir 
is restored no matter how small the circuit radius. Even 
though the degeneracy was lifted by spin-orbit interac- 
tions, the hidden conical intersection exposed by a Zee- 
man field controls the result. 
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